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Abstract 



Exotic dark matter together with the vacuum energy (cosmological con- 
stant) seem to dominate in the flat Universe. Thus direct dark matter de- 
tection is central to particle physics and cosmology. Supersymmetry pro- 
vides a natural dark matter candidate, the lightest supersymmetric particle 
(LSP). It is possible to obtain detectable rates, but realystically they are 
expected to be much lower than the present experimental goals. So one 
should exploit two characteristic signatures of the reaction, namely the 
modulation effect and the correlation with the sun's motion in directional 
experiments. In standard non directional experiments the modulation is 
small, less than two per cent. In the case of the directional experiments, 
the main subject of this paper, we find two novel features, which are essen- 
tially independent of the SUSY model employed, namely: 1) The forward- 
backward asymmetry, with respect to the sun's direction of motion, is very 
large and 2) The modulation observed in a plane perpendicular to the sun's 
motion can be higher than 20 per cent and is direction dependent. 

PACS numbers:95.35.+d, 12.60. Jv 
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I. INTRODUCTION 



The combined MAXIMA-1 [1], BOOMERANG [2], DASI [3] and COBE/DMR Cos- 
mic Microwave Background (CMB) observations [4] imply that the Universe is flat [5], 
Q = 1.11 ± 0.07, while the baryonic component is very small Vt b h 2 = 0.032 + ° og 9 . Fur- 
thermore exotic dark matter has become necessary in order to close the Universe. In 
fact about a decade ago the COBE data [4] suggested that CDM (Cold Dark Matter) 
dominates the Universe, VLcdm being at least 60% [6]. Subsequent evidence from two 
different teams, the High-z Supernova Search Team [7] and the Supernova Cosmology 
Project [8] ' [9] changed this view suggesting that the Universe may be dominated by 
the cosmological constant A or dark energy. In other words one roughly finds a baryonic 
component Qb = 0.1 along with the exotic components VLcdm = 0.3 and Qa = 0.6. In a 
more detailed recent ACDM analysis by Primack [10] (see also the analysis by Einasto 
[11]) we find h = 0.72 ± 0.08, tt b = 0.040 ± 0.002, and Vt m = Vt CDM = 0.33 ± 0.035 (from 
cluster baryons etc), 0.4| ± 0.2 (from cluster evolution) and 0.34 ±0.1 (from Lya forest 
P(k)), while VL HDM < 0.05 and fi A = 0.73 ± 0.08. In other words Q m « 3/4 Q A - Since 
the non exotic component cannot exceed 40% of the CDM [12], there is room for the 
exotic WIMP's (Weakly Interacting Massive Particles). In fact the DAMA experiment 
[13] has claimed the observation of one signal in direct detection of a WIMP, which with 
better statistics has subsequently been interpreted as a modulation signal [14]. 

Supersymmetry naturally provides candidates for the dark matter constituents [15]- 
[18]. In the most favored scenario of supersymmetry the LSP can be simply described 
as a Majorana fermion, a linear combination of the neutral components of the gauginos 
and Higgsinos [15]- [22]. The main ingredients of the SUSY input are summarized in 
the next section. The essential features of the LSP-nucleus cross-section are discussed in 
sect. 3. The basic formulas for the event rates, in the case of Gaussian models as well 
as those due caustic rings, are given in sect. 4. The results obtained both for the non 
directional as well as the directional experiments are discussed in sect. 5. Finally in sect. 
6 we present our conclusions. 
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II. THE ESSENTIAL THEORETICAL INGREDIENTS OF DIRECT 

DETECTION. 



It is well known that there exists indirect evidence for a halo of dark matter from 
the observed rotational curves. It is, however, essential to directly detect [15]- [26] such 
matter, since this, among other things, will also unravel the nature of the constituents 
of dark matter. The possibility of detection depends on the nature of such constituents. 
Here we will assume that such a constituent is the LSP. Since this particle is expected 
to be very massive, m x > 30GeV, and extremely non relativistic with average kinetic 
energy T 50KeV(m x /100GeV), it can be directly detected [15]- [26] mainly via the 
recoiling of a nucleus (A,Z) in elastic scattering. The event rate for such a process can 
be computed from the following ingredients: 

1) An effective Lagrangian at the elementary particle (quark) level obtained in the 
framework of supersymmetry as described , e.g., in Refs [22,23]. 

2) A well defined procedure for transforming the amplitude obtained from the previous 
effective Lagrangian from the quark to the nucleon level, i.e. a quark model for the 
nucleon. This step is not trivial since the obtained results depend crucially on the content 
of the nucleon in quarks other than u and d. This is particularly true for the scalar 
couplings, which are proportional to the quark masses [27]— [30] as well as the isoscalar 
axial coupling. 

3) Compute the relevant nuclear matrix elements [31]— [33]. using as reliable as 
possible many body nuclear wave functions. Fortunately in the case of the scalar coupling, 
which is viewed as the most important, the situation is a bit simpler, since then one needs 
only the nuclear form factor. 

Since the obtained rates are very low, one would like to be able to exploit the modu- 
lation of the event rates due to the earth's revolution around the sun [34,35] [36]— [39]. 
In order to accomplish this one adopts a folding procedure, i.e one has to assume some 
velocity distribution [34,40], [37,39], [41]- [43] for the LSP. One also would like to exploit 
the signatures expected to show up in directional experiments, by observing the nucleus 
in a certain direction. Since the sun is moving with relatively high velocity with respect 
to the center of the galaxy one expects strong correlation of such observations with the 
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motion of the sun [19,44]. On top of this one expects in addition to see a more interesting 
pattern of modulation. 

The calculation of this cross section has become pretty standard. One starts with 
representative input in the restricted SUSY parameter space as described in the litera- 
ture [20,22] (see also Arnowitt and Dutta [24]). We have adopted a similar procedure, 
which has previously described [20] and will not be repeated here. In the above procedure 
the most important constraints on the SUSY parameter space come from [24,20]: 

1. ) The LSP relic abundance which must satisfy the cosmo logical constrain: 

0.09 < tt LSP h 2 < 0.22 (1) 

2. ) The Higgs mass bound. The bound is obtained from the recent experiments CDF 
experiment [45] and, especially the LEP2 signal [46], [47] m h = 115^0 9 GeV, in SUSY 
cannot unambiguously attributed to a definite mass eigenstate. Furthermore in the LSP- 
nucleus scattering both physical scalar eigenstates contribute (the surviving pseudoscalar 
does not lead to coherence). The correct prediction for the Higgs mass, however, is 
essential, since it imposes important constraints on the allowed parameter space [48], 
[49]. 

3.) The 6 ^ s 7 limit and the bound on the anomalous magnetic moment of the 
muon, see, e.g., recent work and references therein [25]. 

4.) The need to will ourselves to LSP-nucleon cross sections for the scalar coupling, which 
gives detectable rates 

4 x 10" 7 pb < a^T < 2 x 10~ 5 pb (2) 

This is because above this range the direct observation of LSP should have occurred in 
the experimental searches so far and below this it is unobservable. We should remember, 
however, that the event rate does not depend only on the nucleon cross section, but 
on other parameters also, mainly on the LSP mass and the nucleus used in target. 
The condition on the nucleon cross section imposes the most severe constraints on the 
acceptable parameter space. In particular in our model [20] as well as in other models 
[51,25] it restricts tanf3 to values tanf3 ~ 50. 
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III. THE LSP-NUCLEUS EXPRESSIONS DIFFERENTIAL CROSS SECTION . 



The expressions for this cross-section are well known. We will, however, summarize 
the main ingredients here for the reader's convenience and in order to establish notation. 

To begin with the effective Lagrangian describing the LSP-nucleus cross section can 
be conveniently put in a form familiar from weak interactions: [23] 

Leff = -^{(xi7 A 75Xi)^a + (X1X1) J} (3) 

where 

Jx = N lx {f v + f v r 3 + f Al5 + f\ lb r,)N , J = N{f s + fS)N (4) 

We have neglected the uninteresting pseudoscalar and tensor currents. Note that, 
due to the Majorana nature of the LSP, Xil X Xi — (identically). 

With the above ingredients the differential cross section can be cast in the form 
[19,36,37] 

da — — it^ — 
d°{u,v) = 2{fibu)2 [^s + F\u) + E spin F n (u)} (5) 

^ = (To( mao ) {A [{fs ~ fs ~ a~ } ] " a »* A ( mao } (6) 

yi _ spin j- 

Z-'spin p,x° ^> s Pi- n ' >>spin jo l ~ > \"') \ • ) 

J A 

S{u) = [(^o(0)) 2 |^ + 2^^0(0)^(0)^ + ^i(O)) 2 ] (8) 
Ja F iii u ) Ja f ii\ u ) 

^ v = a P,x° & ( 9 ) 



(/i r (A)//i r (iV)) 2 f y A-2Z VQ 
& 7. fh7~, A ^ ~ To a — ) K— ) i L 



1 2?7 + 1 ( 2u ). 



1 + 4) 2 fv A c (2/^6)2 (1 + 77)2 ( W 2 



(10) 
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a p x o — proton cross-section, % = S, spin, V given by: 

a px° = a ° (fs) 2 ( M ^ ) 2 (scalar) , (the isovector scalar is negligible, mainly since the 
heavy quarks dominate [27]— [30], i.e. = erf) 

<*$ = °* 3 {fl + f\) 2 {^f (spin) , < x0 = a (ft + f v f {^f (vector) 

where is the nucleon mass, n = m x jm^A, and fi r (A) is the LSP-nucleus reduced 

mass, fJ> r (N) is the LSP-nucleon reduced mass and 

a = }-{G F m N f ~ 0.77 x 10~ 38 cm 2 (11) 

Q = Qou , Q = -^—=A.lxlQ A A-^KeV (12) 

where Q is the energy transfer to the nucleus and F(u) is the nuclear form factor and 
Fu(u) is the isovector spin response factor. S(u) is essentially independent of u. It 
depends crucially on the static spin matrix elements and the ratio of the elementary 
isoscalar to isovector amplitudes [33] . 

In the present paper we will consider both the coherent mode as well as the spin 
mode, but we will not focus on the nucleon cross section. For the scalar interaction we 
will use a form factor obtained as discussed in our earlier work. We will also consider 
the spin contribution, which is expected to be more important in the case of the light 
targets. For a discussion of the spin matrix elements we refer the reader to the literature 
( see Divari et al [33,50].) We only mention here that the spin matrix elements are the 
largest and the most accurately calculated for the A = 19 system. The static spin values, 
however, affect the quantity S spin , which affects the expected rate, but it is not the main 
subject of the present work. What is most relevant here is the spin response function 
F n (u). For the light nuclei it was taken from our earlier work [33] and or the 127 1 target 
it was obtained from the calculation of Ressel et al [31]. 

The vector contribution, which, due to the Majorana nature of the LSP, is only a 
relativistic correction and, at present, below the level of the planned experiments, can 
safely be neglected. 
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IV. EXPRESSIONS FOR THE RATES. 



The non-directional event rate is given by: 

dN p(0) m 

R = Rnon-dir = -7T = ~ A Cr(u,V)\v\ (13) 

dt m x ArriN 

Where p(0) = 0.3GeV/cm 3 is the LSP density in our vicinity and m is the detector mass 
The differential non-directional rate can be written as 

dR = dR non _ dir = m da(u,v)\v\ (14) 

m x ArriN 

where da(u, v) was given above. 

The directional differential rate [19], [39] in the direction e is given by : 

O ( ) Tfl \ 

dRdir = - t v.eHv.e) — da(u,v) (15) 

m x Am N Z7r 

where H the Heaviside step function. The factor of l/2n is introduced, since the differen- 
tial cross section of the last equation is the same with that entering the non-directional 
rate, i.e. after an integration over the azimuthal angle around the nuclear momentum has 
been performed. In other words, crudely speaking, l/(27r) is the suppression factor we 
expect in the directional rate compared to the usual one. The precise suppression factor 
depends, of course, on the direction of observation. The mean value of the non-directional 
event rate of Eq. (14), is obtained by convoluting the above expressions with the LSP 
velocity distribution f(v,VE) with respect to the Earth, which moves with velocity ve 
relative to the sun (see below), i.e. is given by: 

/dR\ p(0) m f da(u,v) 3 

The above expression can be more conveniently written as 



Now we perform the needed integrations. First over the velocity distribution ranging 
from avo^/u, where a = \\/2mu r bvQb]~ l and b the nuclear (harmonic oscillator) length 
parameter, to the maximum escape velocity v m . Second over the energy transfer u 

6 



ranging from u min dictated by the detector energy cutoff to u max = (v m /(v a)) 2 . Thus 
we get: 



where t is the total rate in the absence of modulation, a is the phase of the Earth 
(a — around June 2nd) and Q m in is the energy transfer cutoff imposed by the detector. 
In the above expressions R is the rate obtained in the conventional approach [23] by 
neglecting the folding with the LSP velocity and the momentum transfer dependence of 
the differential cross section, i.e. by 



where Sj, i — S, V, spin contain all the parameters of the SUSY models. The modulation 
is described by the parameter h . 

The total directional event rates can be obtained in a similar fashion by suitably 
modifying Eq. (16) 



The role played by the velocity distribution is very clear. What is the proper velocity 
distribution to use? The best approach seems to be to apply the Eddington approach [53]. 
In such an approach one starts from a given density as a function of space and one solves 
Poisson's equation to obtain the potential. Then from the functional dependence of the 
density on the potential one can construct, at least numerically, the density distribution 
in phase space, as a function of the potential $(r) and the velocity [53], [41,42,54]. 
Evaluating this distribution in our vicinity yields the desired velocity distribution. Since 
this procedure can only be implemented numerically it is very hard to incorporate in 
the calculation of the directional rates. We thus follow the conventional approach and 
use two phenomenological velocity distributions: i) A Gaussian distribution, which can 
be either spherically symmetric or only axially symmetric and ii) A velocity distribution 
prescribed by the assumption of the late infall in the form of caustic rings. 



R — R t [1 + h(a, Q min )cosa)] 



(18) 




(19) 




(20) 
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A. Gaussian Distribution 



The Gaussian distribution, with respect to the center of the galaxy, is of the form: 

(A+1)[K) 2 + K) 2 ] + K) 



f(v,\)=N(\) 7 -^=-Exp 



(21) 



where vq is equal to the velocity of the sun around the center of galaxy, A is the asymmetry 
parameter, assumed to be in the range < A < 1, and N(X) is a normalization constant, 
iV(0) = 1. One must of, course, transform the above distribution into the local coordinate 
system, taking into account both the motion of the sun as well as that of the Earth. In 
writing the above velocity distribution we have chosen a set of axes as follows: 

The z-axis is along the sun's direction of motion. 

The x-axis is the radial direction outwards. 

The y- axis is perpendicular to the galactic plane, so that the system is a right-hand 
one. 
Then 

v =v + v z + ve, ve = ve [sinarr — cos a cos 7?/ + cos a sin 7^] (22) 

with 7 7r/6 and a the phase of the Earth. After that one must do the folding with the 
above velocity distribution. The integration, however, of equation 20 is quite difficult 
due to the presence of the Heavyside function. So for the purpose of integration we found 
it convenient to go to a coordinate system in which the polar axis is in the direction of 
of observation e, which in the above coordinate system is specified by the polar angle 
and the azymouthal angle $. In this new coordinate system polar angle specifying 
the velocity vector is simply restricted to be < 9 < n, while the azymouthal angle 
is unrestricted. Thus the unit vectors along the new coordinate axes, X, Y, Z, are 
expressed in terms of the old ones as follows: 

Z = sin cos + sin sin $y + cos Qz. (23) 
X = cos cos $x + cos sin $y — sin Ox. (24) 
Y = - sin + cos (25) 
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Thus the LSP velocity is expressed in the new coordinate system as: 

v x = sin 6 cos §vx + sin 6 sin Qvy + cos Qvz, (26) 

v y = cos 6 cos &v x + cos 6 sin §v Y — sin <8>v z , (27) 

v z = - sin §v x + cos &v Y , (28) 

with vx = fsin#cos0, -uy = -u sin 6* sin 0, vz = vcosO. It is thus straightforward to go 
to polar coordinates in velocity space and get: 

\ du l dir m x Amu Jav ^i Jo Jo 2tt du 

with ^ = cosO. Now the orientation parameters 9 and $ appear explicitly in the distri- 
bution function and not implicitly via the limits of integration. The function / can be 
obtained from the velocity distribution, but it will not be explicitly shown here. Thus 
we obtain: 

Rdir = R(tdir/2n) [1 + h m cos(a - a m n)} (30) 

where the quantity tdi r provides the non modulated amplitude, while h m describes the 
modulation amplitude and a m is the shift in phase (in units of n), giving the phase of the 
Earth in which the maximum modulation occurs. Clearly these parameters are functions 
of and $ as well as the parameters a and Q m in- The dependence on a comes from the 
nuclear form factor the folding with LSP velocity. The other SUSY parameters have all 
been absorbed in R. 

Instead of tdir itself it is more convenient to present the reduction factor of the non 
modulated directional rate compared to the usual non-directional one, i.e. 

fred = ^jr = W(2?r t) = k/(2tt) (31) 

It turns out that the parameter k, being the ratio of two rates, is less dependent on the 
parameters of the theory. The directional rate can be cast in an even better form as 
follows: 
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R dir = R t— [1 + h m cos(a - a m tt)] (32) 

Another quantity, which may be of experimental interest is the asymmetry As = 
\Rdir{— ) — Rdir(+)\/ (Rdir(-) + Rdir(+)) in some given direction (+) and its opposite 
(-). The most relevant direction for the asymmetry is that of the velocity of the sun. 
As is almost independent of all other parameters except the direction of observation and 
the velocity distribution. The directional rates exhibit interesting pattern of modulation. 
From the functions h m (a, Q m in) and a m obtained in the present work, it is trivial to plot 
the the expression (32) as a function of the phase of the earth a. Hence, this will not be 
done here. 



B. Caustic Rings 

One would like to examine a non isothermal model to see what effect, if any, may 
have on the directional rate. The model of caustic rings proposed by Sikivie [55] comes 
to mind [56], [39]. Admittedly, however, this scenario for dark matter distribution in the 
galaxy is not broadly accepted [57], since it has not been supported by galaxy evaluation 
simulations. . From previous work [39] we take the data for caustic rings needed for our 
purpose , which, for the reader's convenience, we have included in Table I. 

The above expressions simplify significantly in the case of caustic rings whereby the 
velocity distribution is discreet. Thus Eq. (17) takes the form 

<^> = ^fy a 2 \E s F (u) + i^ixvAH + S spm F spm ( M )] (33) 

We remind the reader that p was obtain for each type of flow (+ or -), which explains 
the factor of two. in the Sikivie model [55] we have (2p/p(0) = 1.25. The quantities 
F ,F 1 , F spin are obtained from the corresponding form factors via the equations 

F k (u) = F 2 {u)^ k {u)^^ , A; = 0,1 (34) 

Fspin 

(u) = F n («)* («) (35) 

with 
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/9 N i, 2 
**(«) = E PnVltM^-u) (36) 
V n=i a 



with 



Van = [(Vnz ~ 1 ~ ^ C ° Sa ^ + ^ n V + ^ C ° S7 C0Sa ) 2 + (Vnx ~ ^ Sma) 2 ] 1/2 (37) 

with 5 = 2{v E /wpsilono). Integrating numerically Eq. (33) we obtain the total undi- 
rectional rate as a function of the phase of the Earth. Unlike our previous work we did 
not make an expansion in terms of 5, in order to better deal with threshold effects. By 
making a Fourier decomposition of the obtained rate to an accuracy better than 1% we 
find : 

R = Rt-^-[l + h(a,Q min )cosa)} (38) 
P(0) 

In other words the modulation is again described a parameter h, which, of course, takes 
different numerical values compared to those of the Gaussian distribution. Note that, 
unlike our earlier work [39], the modulation h in this work is defined with a + sign. So 
our present value of h in the case of caustic rings is expected to have a phase difference 
of 7r compared to that of the usual Gaussian distribution. 
Similarly integrating Eq. (20) we get 

dT, 2o a 2 (v 2 ) - 

(^)dir = —[Z s F (u) + A_/E yFl ( u ) + E spin F spin (u)} (39) 
The quantities F , F 1? F spin are obtained from the equations 

F k (u) = F\u)^ k (u)^±^-, k = 0, 1 (40) 

F spin (u) = Fn(«)*o(«) (41) 

N 

**(«)= rE PnvT^XHiX) (42) 

71=1 

with H(X) the usual Heaviside (theta) function and X given by 

5 5 5 

X = (y nz - 1 - - siwy cosa)e z .e + (y ny + - cosj cosa)e y .e + (y nx - - sina)e x .e (43) 

11 



Note that in X only certain rings contribute for a given direction of observation e (as 
dictated by the Heaviside function). Note further that the other Heaviside function for 
a given ring n restricts the contribution of the form factor as follows: 



and an analogous expression for the spin response function. 

Once again we did not make an expansion in powers of delta. We Fourier decomposed 
the final expression and we found that Fourier components higher than unity (n > 2) are 
negligible. In other words our results for the directional rate can be cast in the form of 



The three basic ingredients for the event rate for the LSP-nucleus scattering are the 
input SUSY parameters, a quark model for the nucleon and the velocity distribution 
combined with the structure of the nuclei involved. In the present work we will present 
our results for the coherent scattering and make comparisons with the spin contribution, 
whenever possible. We will focus our discussion on the light targets (e.g. A=19, 23, 29) 
and the more popular target 127 1. 

We have utilized two nucleon models indicated by B and C, for their description see 
our previous work [50], which take into account the presence of heavy quarks in the 
nucleon. We also considered the effects on the rates of the energy cut off imposed by the 
detector, by considering two typical cases Q m % n = 0, 10 KeV . 

The experimentally determined quenching factors, see Simon et al and Graichen et 
al [58,59] and references therein, defined as the ratio of a signal induced by nuclear 
recoil to that of an electron of the same energy, have not been included in calculating 
the total rates in the present work. These factors, which are functions of the energy, 
depend mainly on the detector material. For Nal they have been measured [58,60] down 
to about 10 KeV and they are constant, about 0.25. For our calculations employing 
Qmin = 10 KeV this amounts to a reduction of the parameter t to about 25% of its 
value presented here. We cannot estimate what the effect of quenching is going below 




(44) 



Eq. (32). 



V. RESULTS AND DISCUSSION 
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10 KeV. The modulation amplitude and the reduction factor k, however, being relative 
parameters, are not expected to be influenced very much by such effects. 

A. Isothermal Models 

For the reader's convenience we will begin by presenting results for the unmodulated 
non directional event rates, Rt, of the symmetric isothermal model for a favorable SUSY 
parameter choice with large tan((3), described previously [20], shown in Fig. 1 as a 
function of the LSP mass in GeV. Clearly, depending on the SUSY parameter space and 
in particular the LSP mass, the rates can change many orders of magnitude. What we 
rally want to exhibit is the role of the nucleon model employed and the effect of the 
energy cutoff. 

The two relative parameters, i.e. the quantities t and h, for 127 1 are shown in Figs 
2 and 3 respectively in the case of Gaussian distribution. For the light systems these 
quantities are essentially constant independent of the LSP mass(the reduced mass does 
not change as the LSP mass increases). They are shown in Table II. In the case of the 
directional rates we calculated the reduction factors and the asymmetry parameters as 
well as the modulation amplitude as functions of the direction of observation, focusing 
our attention along the three axes [38], i.e along +z, —z, +y, —y, +x and —x. In the case 
of the directional rates we calculated the reduction factors and the asymmetry parameters 
as well as the modulation amplitude as functions of the direction of observation, focusing 
our attention along the three axes [38], i.e along +z, —z, +y, —y, +x and —x. 

Since f re d is the ratio of two parameters, its dependence on Q m in and the LSP mass is 
mild. In the case of light nuclear systems (A=19, 23 and 29) these parameters are shown 
in Table II. Note that in the favorable direction —z (opposite to the velocity of the sun) 
the modulation is about a factor of two bigger than in the non-directional case, but it 
is still quite small (h = 0.06). The reduction factor is k — 0.7. In the sun's direction of 
motion the process is unobservable. In the plane perpendicular in the sun's velocity the 
rate is further reduced by a factor of about three. 
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FIGURES 
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FIG. 1. The Total detection rate per (kg — target)yr vs the LSP mass in GeV for a typical 
solution of the parameter space, as described in our previous work (see text), in the case of 
127 1. The results shown by thick lines correspond correspond to model B, while those by a thin 
line to model C. In the upper curve no detector cutoff was employed, while in the lower curve 
we used a detector energy cutoff of Q m in = 10 KeV. Such effects introduce variations in the 
rates by factors of about two. 
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FIG. 2. The dependence of the quantity t on the LSP mass for the symmetric case (A = 0) 
on the left as well as for the maximum axial asymmetry (A = 1) on the right in the case of 
the intermediate mass target 127 1. For orientation purposes two detection cutoff energies are 
exhibited, Q m in = (thick solid line) and Q m in = 10 keV (thin solid line). As expected t 
decreases as the cutoff energy and/or the LSP mass increase. We see that the parameter A has 
little effect on the non modulated rate. 
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FIG. 3. The same as in Fig. 2 for the modulated amplitude. We see that the modulation 
is small and decreases with the LSP mass. It even changes sign for large LSP mass. The 
introduction of a cutoff Q rn in increases the modulation (at the expense of the total number of 
counts). It also increases with the parameter A. 

The modulation shows a very interesting pattern. If the observation is done in the 
direction opposite to the sun's direction of motion, the modulation amplitude h m behaves 
in the same way as the non directional one, namely h. It is more instructive to consider 
directions of observation in the plane perpendicular to the sun's direction of motion 
(0 = tt/2). We see from Table II that now the modulation is quite a bit larger, giving 
rise to a difference between the maximum and the minimum rates of about 50%. 

For heavier nuclei the pattern slightly changes. Our results are presented in plots, 
in which we adopted the following convention: The thick solid line corresponds to +z 
(0 = 0), while the thin line to —z (0 = tt. In the case of = tt/2 we encounter 4 cases. 
The intermediate thickness line corresponds to the ±x, the dotted line to +y and the 
dashed line to —y. In some cases two or more lines may coincide. In the case of k one 
can distinguish only the curves corresponding to the three values. 
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The quantities k and h m show some variation with the LSP mass (see Figs 4 and 6), 
since in this case the reduced mass changes. 

The quantities As and a m do not show any significant changes compared to those 
of the light systems (see Fig. 5 and 7. We see that, in the absence of modulation, the 
asymmetry is non zero only if e is in the direction of the sun's motion. In the other 
directions the asymmetry depends on the phase of the Earth and is the same with the 
modulation , i.e. h m \ cosa| , h m \ sina\ in the y-direction (perpendicular in the plane of 
the galaxy) and x-direction (in the radial direction in the galaxy) respectively. 

k hm 
1 0.5 

0.8 0.4 




. 2 - . 1 

25 50 75 100 125 150 175 2o8' 1SP ~ > 25 50 75 100 125 150 175 2oS llSP 



FIG. 4. On the left figure one sees the quantity k and on the right the quantity h m for 
A = and Q m in = 0. The results are almost identical for the coherent and the spin modes. 
They change very little for Q m i n = 10 KeV. For the identification of the curves see text. Note 
that in the case of the modulation the curves corresponding to +z and ±y coincide. The large 
modulation seen in the +z direction is essentially useless since the event rate is very small. 
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0.6 
0.4 
0.2 



am 
1 . 6r_ 

1.4 

1.2 



25 50 75 100 125 150 175 20 



0.2 



nlsp-> 



25 50 75 100 125 150 175 20 



nlsp-> 



FIG. 5. The quantities As and a m for A = and Q m in = 0. The asymmetry As, shown 
on the left, takes only two values As ~ 1.0 in the direction of the sun's motion and zero in all 
other directions. These asymmetry plots do not contain the contribution due to modulation, 
since, then, the asymmetry would depend on the time of observation (see text). On the right 
we show the shift (in units of ir) in the position of the maximum of the modulated amplitude. 
This shift is almost zero in the — z, +y directions, close to ir in the +z, —y directions, close to 
7r/2 in the +x direction and almost 3n/2 in the — x direction. The notation is the same as in 
Fig. (4). 
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0.4 
0.2 
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25 50 75 100 125 150 175 20 
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FIG. 6. The same as in Fig. 4 for A = 1. 
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am 

1 . 6r 

1 . 4 
1 . 2 



0.8 
0.6 
0.4" 
0.2 

50 75 100 125 150 175 200 1SP ~ > 25 50 75 100 125 150 175 200 1SP ~ > 

FIG. 7. The same as in Fig. 5 for A = 1. 
The results for the light systems are presented in Table II. 

B. Caustic Rings 

The case of non isothermal models, e.g. caustic rings, has previously been discussed 
[39] in some special cases. Here we will expand our discussion further. In the case of the 
non directional rates, the quantities t and h for 127 1 are shown in Figs 8. We see that as 
far as t is concerned there is no essential difference between caustic rings and Gaussian 
distribution. Notice, however, that the modulation is now smaller and of opposite sign. 
In the case of directional rates the quantities k and h rn for 127 1 are shown in Figs 9 and 
10 for Qmin = 0, 10 KeV respectively. The results obtained for the light systems are 
shown in Table III. One clearly sees that the maximum rate is now in the direction 
of the sun's motion, +z, and the minimum in the — z direction, i.e. in the opposite 
sense compared with the Gaussian distribution. In the other directions the rates fall in 
between. Naturally the rate is reduced in the presence of an the energy cutoff, but in the 
case of caustic rings the reduction manifests itself mainly for small LSP masses. For such 
masses it is not easy to have energy above threshold transferred to the nucleus, motion, 
but smaller than in the Gaussian distribution (As = 0.75 and As = 0.68 for Q mm = 
and 10 KeV respectively). In the other directions the asymmetry is governed by the 
modulation. 



As 
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FIG. 8. The same as in Fig. 2 (left figure) and 3 (right figure) in the case of caustic rings. 
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FIG. 9. The same as in Fig. 4 for caustic rings with Q m %n = (A is irrelevant). 




19 



VI. CONCLUSIONS 



In the present paper we have discussed the parameters, which describe the event 
rates for direct detection of SUSY dark matter. Only in a small segment of the allowed 
parameter space the rates are above the present experimental goals [20,22,24]. We 
thus looked for characteristic experimental signatures for background reduction, i.e. a) 
Correlation of the event rates with the motion of the Earth (modulation effect) and b) 
the directional rates (their correlation both with the velocity of the sun and that of the 
Earth.) 

A typical graph for the total non modulated rate is shown Fig. 1. The relative 
parameters t and h in the case of non directional experiments are exhibited in Fig. 2 and 
Fig. 3 for the Gaussian models. For caustic rings they are shown in Fig. 8. We must 
emphasize that the two graphs of Figs. 2 and 8 do not contain the entire dependence 
on the LSP mass. This is due to the fact that there is the extra factor m^ 1 in Eq. (19) 
and a factor of $ arising from Sj, i = S, spin, V (see Eqs (6), (7), and (10)). All these 
factors combined lead to a essentially a constant. There remains, however, an LSP mass 
dependence, which is due to the fact that the nucleon cross section itself dramatically 
depends on the LSP mass. 

Figs 2, 3 and 8 were obtained for the scalar interaction, but similar behavior is 
expected for the spin contribution. From the point of view of the static spin matrix 
elements the most favored system is the A = 19 (see our previous work [33]). The 
scale of the total spin contribution, coming from the SUSY dependent parameter E s , 
which was not discussed in this work, may be very very different from the corresponding 
to the scalar amplitude quantity. We should also mention that in the non directional 
experiments the modulation 2h is small, .i.e. for A = less than 4% for Qmin = and 
increases to 12% for Qmin = 10 KeV (at the expense of the total number of counts). For 
A = 1 there in no change for Q m in = 0, but it can go as high as 24% for Q m in = 10 KeV. 
In the case of caustic rings is smaller and of opposite sign. 

For the directional rates it is instructive to examine the reduction factors k if 
the observation is made in a specific direction, e.g. along the three axes, i.e along 
+z, —z, +y, —y, +x and —x. These depend on the nuclear parameters, the reduced mass, 
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the energy cutoff Q m in an d A [38]. Since f re d is the ratio of two parameters, its dependence 
on Qmin and the LSP mass is mild. So we present results for 127 1 in Figs 4 and 6 (see also 
Figs 9 and 10 for caustic rings). In the case of light systems our results are presented in 
Tables II and III. As expected the maximum rate is along the sun's direction of motion, 
i.e opposite to its velocity (—z) in the Gaussian distribution and +z in the case of caustic 
rings. In fact we find that k(—z) is around 0.5 (A = 0) and around 0.6 (A = 1.0). It 
is not very different from the naively expected f re d = l/(27r) , k — 1. The asymmetry 
along the sun's direction of motion, As = \Rdir(— ) — Rdir(+)\/(Rdir(—) + Rdir(+)) is 
quite characteristic, i.e. almost unity for Gaussian models and a bit smaller in the case 
of caustic rings. The rate in the other directions is quite a bit smaller (see Tables II and 
III) and the asymmetry is equal to the absolute value of the modulation. 

The disadvantage of smaller rates in the plane perpendicular to the sun's velocity 
may be compensated by the bonus of very large and characteristic modulation. 

In conclusion: in the case of directional non modulated rates we expect unambiguous 
correlation with the motion on the sun, which can be explored by the experimentalists. 
If one concentrates in a given direction there appears a reduction factor. The reduc- 
tion factor in the most favored direction, .i.e in the line of the motion of the sun, is 
approximately only 1 / (4tt) relative to the non directional experiments. In the plane per- 
pendicular to the motion of the sun we expect interesting modulation signals, but the 
reduction factor becomes worse. These reduction factors do not appear to be an obstacle 
to the experiments, since the TPC counters to be used in the planned experiments can 
make observations in almost all directions simultaneously [61]. Once some candidate 
events are seen, one can analyze them further in the way we propose here by selecting 
those corresponding to a given direction of observation. 

This work was supported by the European Union under the contracts RTN No HPRN- 
CT-2000-00148 and TMR No. ERBFMRX-CT96-0090. During the last stages of this 
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TABLES 

TABLE I. The quantities a n ,yi' n = v n /v ,y nz = v n< f,/v ,y ny = v nz /v ,y nx = v nr /v and 
Pn = d n /[J2^Li d n ] and y n = [(y nz - l) 2 + y 1 + yl x ] 1/2 (for the other definitions see text ). 
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TABLE II. The quantities K,h rn ,a m and As for the light nuclear systems (A=19, 23 and 
29), are almost identical. The difference between the coherent and the spin mode is less than 
1%. They are also essentially independent of the LSP mass and the SUSY parameters. Note 
the phase of the modulation in the directions —z, + y is the same as in the non directional 
case. The phase of the modulation in the — y direction is reversed (minimum in June 3nd). The 
phase in the ±x directions leads to a maximum in between. When the range of a variable is 
given, it depends somewhat on the LSP mass (the mass increases to the right). 
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quantity dir 



A = 



A = A = 1.0 



A = 1.0 







Qmin — 0.0 


Q min = 10.0 KeV 


Qmin — 0.0 


Qmin = 10.0 KeV. 


t 


all 


1.0 


0.6 


1.0 


0.6 


h 


all 


0.02 


0.04-0.06 


0.04 


0.05-0.08 




+z 


0.018 


0.010 


0.003 


0.000 




X 


0.190 


0.177 


0.211 


0.145 - 0.186 


K 


y 


0.190 


0.177 - 0.180 


0.150 


0.087 - 0.125 




-z 


0.690 


0.758 - 0.760 


0.752 


0.519 - 0.672 




+z 


0.211 


0.242 - 0.226 


0.361 


0.380 




+x 


0.235 


0.292 - 0.255 


0.237 


0.325 - 0.261 




+y 


0.199 


0.299 - 0.233 


0.290 


0.456 - 0.347 


h m 


-X 


0.235 


0.292 - 0.255 


0.237 


0.325 - 0.261 




-y 


0.199 


0.199 - 0.208 


0.290 


0.243 - 0.280 




-z 


0.060 


0.100 - 0.068 
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0.158 - 0.092 




+z 
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1 


1 


1 




+x 


1/2 


0.445 - 0.484 


1/2 


0.432 - 0.467 




+y 
















-X 


3/2 
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3/2 


1.587 - 1.533 




-y 


1 


1 


1 


1 




-z 
















z 


0.945 


0.989 - 0.970 


0.991 


1.000 


As 


X 


h m \ sin a 


h m \ sin a 


h m \ sin a 


h m \ sin a 


As 


y 


h m \ cos a 


h m \ cos a 


h m \ cos a 


h m \ cos a 
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TABLE III. The same as in Table II in the case of caustic rings for for the light nuclear 
systems (A=19, 23 and 29). Note the difference in the phases between the modulations of the 
two tables. 
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quantity dir coherent 



coherent 



spin 



spin 







Qmin — 0.0 


Q min = 10.0 KeV 


Qmin — 0.0 


Qmin = 10.0 KeV. 


t 


all 


1.23 - 1.20 


0.547 - 0.893 


1.29 - 1.14 


0.520 - 0.845 


h 


all 


-0.015 


-0.057 - (-0.026) 


- 0.015 


-0.057 - (-0.025) 




+z 


0.381 


0.310 - 0.364 


0.383 


0.310 - 0.363 




X 


0.297 


0.281 - 0.293 


0.297 


0.281 - 0.293 


K 


y 


0.150 


0.182 - 0.159 


0.149 


0.183 - 0.159 




-z 


0.060 


0.067 - 0.062 


0.060 


0.068 - 0.062 




+z 


0.086 


0.208 - 0.112 


0.083 


0.310 - 0.113 




+x 


0.143 


0.351 - 0.195 


0.139 


0.351 - 0.195 




+y 


0.249 


0.239 - 0.246 


0.249 


0.239 - 0.246 


h m 


-X 


0.144 


0.351 - 0.195 


0.137 


0.351 - 0.195 




-y 


0.247 


0.315 - 0.267 


0.245 


0.281 - 0.267 




-z 


0.178 


0.209 - 0.187 


0.178 


0.209 - 0.187 




+z 


1 


1 


1 


1 




+x 


3/2 


3/2 


3/2 


3/2 




+y 
















-X 


1/2 


1/2 


1/2 


1/2 




-y 


1 


1 


1 


1 




-z 
















z 


0.726 


0.667 


0.728 


0.667 


As 


X 


h m \ sin a 


h m \ sin a 


h m \ sin a 


h m \ sin a 


As 


y 


h m \ cos a 


h m \ cos a 


h m \ cos a 


h m \ cos a 
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amHmlspL 
1.6 p. 

1 . 4 

1.2 

1 

. 8 
0.6 
0.4 
0.2 



mlsp-> 

50 75 100 125 150 175 200 ^ 



